Superradiance from crystals of molecular nanomagnets 
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We show that crystals of molecular nanomagnets can exhibit giant magnetic relaxation due to the 
Dicke superradiance of electromagnetic waves. Rigorous treatment of the superradiance induced by 
a field pulse is presented. 

PACS numbers: 75.50.Xx, 42.50.Fx 
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High-spin molecular nanomagnets, such as spin-10 
Mni2 and Fes, represent the boundary between classi- 
cal and quantum physics. On one hand, they exhibit 
pronounced magnetic hysteresis, as classical magnets do. 
On the other hand, the very same magnetization curve re- 
veals quantum nature of spin jjj . These unique properties 
of molecular nanomagnets are consequence of long-living 
metastable spin states || due to the large value of spin 
and high energy barriers. The lifetime of these states is 
believed to be dominated by spin-phonon processes [|[ ||] 
and by quantum spin tunneling J|, ^, [?], |8|. The latter, 
for a field-sweep experiment, has been successfully de- 
scribed in terms of single-molecule Landau-Zener (LZ) 
transitions §, (fi], |l| |l| ^ ||[ @. 

Recent ESR experiments 0, |l|, |l|, f|], || have 
demonstrated noticeable resonant absorption of electro- 
magnetic radiation by molecular magnets. In this Letter 
we show that crystals of magnetic molecules can also be 
a powerful source of coherent electromagnetic radiation. 
At low fields, the magnetic relaxation of molecular nano- 
magnets may be affected by the distribution of energy 
levels due to dipolar fields j2^] , nuclear spins |2^, , and 
crystal defects |l3| . Here we will study the case when the 
resonant tunneling of the spin of individual molecules is 
dominated by the large external magnetic field (or by 
a large transverse anisitropy). We will be interested in 
relaxation between the tunnel-splitted ground states of 
magnetic molecules on the two sides of the barrier which 
has two specific features. Firstly, the change of the pro- 
jection of the spin on the symmetry axis, S z , is of order 
S > 1 and the coupling with the light is increased by 
a factor of S in comparison with the usual processes in 
which S z changes by one. Secondly, phonon relaxation 
processes between these states appear only at high orders 
of the perturbation theory and are strongly suppressed 
[B4j hence electromagnetic relaxation will prevail. 

One might think that the individual relaxation of 
molecules is always a good approximation. Note, how- 
ever, that in molecular nanomagnets the wavelength A of 
the electromagnetic radiation due to transitions between 
spin levels is typically greater than the dimensions of the 
crystal L. Consequently, the molecules can coherently 
interact through the electromagnetic radiation that they 
are emitting, thus greatly enhancing its intensity |2q| . 
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FIG. 1: A pair of tunnel-splitted levels vs. energy bias 
W. Coherent light is emitted after crossing the resonance via 
superradiance. 



Consider a crystal of molecular magnets (e.g., Mn 12 or 
Feg) of spin S, with a large energy barrier between spin- 
up and spin-down states, near a resonant value of the 
magnetic field H z . Since the tunnel splitting A is small 
compared to other energy distances, the pair of splitted 
energy eigenvalues can be approximated by (see Fig. [I]). 



: = - ( c, +ei± VW 2 + A 2 



(1) 



where W = e% — 62, and the eigenstates are given by 
1 



(2) 



with C± = J 1 ± IF/VA 2 + W 2 Q. Here Vi and vj 2 
are eigenstates in the left and right wells with tunneling 
neglected. If the transverse field H x is small, the states 
and ip2 are just \S) and — \S). In the strong field the 
ground states of the system in the wells are described by 
the tilted spin and its projection (S z )q. In this case 



W = 2(S z ) gfi B H z 



(3) 



For S ^> 1 the classical approximation for (S z )q works 
well. In particular, for Mni2 (Ti = —DS Z — gfiBH x S x ) 
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FIG. 2: Energy levels of the Mni2 Ac molecule for different 
transverse fields H x off resonance (qhbHz = D/2), labeled 
by (Sz). States with (S z ) ~ are not localized in one of the 
potential wells. h z ~ 0.63 (H x ~ 6T) generates the ground- 
state splitting A/ks 2± 0.15K. 

one has 

(S z ) « Sy/l-h*, h x = g^ B H x /(2SD). (4) 

To create an appreciable tunnel splitting, one has to ap- 
ply a strong transverse field (e.g., H x ~ 6T for A/k B = 
0.15 K for Mni2). For such values of H x potential wells 
become shallow (see Fig. ||) and most of the levels are no 
longer localized in the wells; still (S z )o w 0.777S' is large. 

We assume that the crystal is in thermal equilibrium 
with the electromagnetic radiation and at first compute 
the rate of transition between ip- and ip+ for a single 
molecule. Interaction with the quantized radiation field 
HpM is V = — g^sHpht.S (Notice that only the 

z-component of H p ht couples to the tunneling system.) 
The decay rate is given by the Fermi golden rule 

Ti = ^|(^ + |^|^)| 2 (5MB) 2 fc 3 coth (Jj?L?) , (5) 

where k — nuj Jc is the wave vector of a photon, n is the 
refraction index of the crystal, and 

w = Va 2 + W 2 /h (6) 
is the photon frequency. The matrix element squared is 

\& + \S z \rl>-)\ 2 = (S z )l A2 ^ w2 , (7) 

making Ti proportional to wA 2 . When k B T < A and 
A/h is as high as 2xl0 10 s" 1 (i.e., A/k B =- 0.15 K), 
Eq. (|) at W < A gives Ti - 7.4 x lO^^s" 1 , which is 
negligible compared with typical experimental transition 
rates. The situation changes drastically, however, when 
one studies interaction of a macroscopic number N of 
magnetic molecules with the electromagnetic radiation. 

When the wave length of emitted photons exceeds the 
linear size of the sample, X^S> L, then the phase of these 



photons is the same throughout the sample and the pro- 
cess of radiation becomes coherent. Near the tunnelling 
resonance the system can be described by the effective 
Hamiltonian 

N 

T-Lcs = -^flHcff • R + Vcb, R = ^t 4 (8) 

i=l 

where Tj is a fictitious spin 1/2 describing the resonant 
states of the z-th molecule, 

3At B H cff = Ae^ + We z , (9) 

V c b = 2(S z )og^BRzH p ht, z is the spin-photon interaction. 
It can be checked that TL c s conserves quantum number 
R 2 = R(R+ 1). We are interested here in the states with 
the maximal value of R, R — i? ma x = N/2. Such states 
are produced when we begin with all magnetic molecules 
in the left well in Fig.l, that is, magnetized down, and 
then sweep the applied field across the resonance. The 
dynamics of the system is that of a large spin R coupled 
to photons, and it is convenient to label its states by the 
quantum number M which is the projection of R on H e a 
of Eq. (||) . The matrix elements of the Hamiltonian Eq. 
(||) for transitions between \R, M) and \R, M± 1) con- 
tain factors ,J(R± M)(R^ M + 1) that are additional 
to those for a single molecule. The latter makes the spon- 
taneous radiation power / to scale as (R+M)(R— M+l). 
If \M\ < R and i? ~ N, then I oz N 2 , which is N times 
the incoherent radiation from a system of N molecules. 
This is the superradiance discovered by Dicke |25j . 

At finite temperatures, the time evolution of a sys- 
tem of tunneling molecules in equilibrium with the ra- 
diation field is described by the density-matrix equation 
(DME) for a spin R. For R 3> 1 and low temperatures, 
Rgfj,BH c g/ '(ksT) = RTiu / '(ksT) S> 1 (but not necessar- 
ily huj/(k B T) > 1, see Ref. Q), this DME goes over to 
the classical Landau-Lifshitz equation, 

h = 7[n x H c ff] — «7[nx [nxH cff ]], n = R/i?, (10) 

where 7 = gun/h is the gyromagnetic ratio, 

is dimcnsionlcss damping coefficient, Ti is given by Eq. 
(|) and H off is defined by Eq. (H). 

The first term in Eq. ( |io[ ) gives dissipationless Landau- 
Zener transitions when the field is swept through the res- 
onance such that W = W(t) satisfies W(±oo) = ±00, 
and the initial condition is n(— 00) = — e z . Indeed, the 
Schrodinger equation for a two-level system is equivalent 
to the equation for a precessing spin. The probability 
Pit) for a molecule to stay in the initial state is given by 

P(t) = [l-n z (t)}/2. (12) 

For W(t) = vt, one obtains the LZ result §, 

P(oo) = P = exp[-7rA 2 /(27iw)]. (13) 



3 



The second term in Eq. (10) describes Dicke superra- 
diance since a contains the large factor R — N/2. Eqs. 
(0), (§, (0) and relation gfJ, B H eB = VA 2 + W 2 yield 
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a = l -N{S z ) 2 g 2 n 3 (t- 
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(14) 



where e is electron charge, m e is electron mass and c is 
the speed of light. Note that a is independent of the 
energy bias W. The corresponding longitudinal decay 
rate is 



T = 2ajH eS = 2aoj = 2RT" 



(o) _ NT (o) 



(15) 



where is given by Eq. (^|) without the coth factor. 
Thus superradiance boosts the one-molecule relaxation 
rate by a macroscopically large factor N. 



Let us consider at first the solution of Eq. (10) for 
H e ff = const. In this case it is convenient to describe the 
motion of the Dicke pseudospin n in the coordinate sys- 
tem with the z' axis along H c ff. The physically relevant 
z component of n is given by 



n z = n z i cos-0 — n x i sin?/>, cos^ 



W 



, (16) 



VA 2 + W 2 

whereas the solution for n z > and n x > reads 

sinh(r) + n z < (0) cosh(r) 
n z ,{t) = )\ )' , r ee a 7 H eS t, 

cosh(r) + n z > (0) smh(r) 

n x ,{t) = ^Jl-n 2 z ,(t) sin(7 H cS t + cf> ). (17) 

The full relaxation from n z i = — 1 to n z i = 1 is given by 
n z i (t) = tanh(r) which can be used to obtain the value 
of a from an independent macroscopic argument. To this 
end, we set B = N/2 and equate the change of Zceman 
energy iVg/isHeff [see Eq. (||)] to the energy dissipated 
due to the magnetic-dipole radiation I(t)dt, 

I=[2/(3c 3 )]m 2 z , (18) 

where m z = N (S z )ogfiBn z . For a <§; 1 one obtains 

n z = -h x > sin^ = {"/H cS ) 2 n x > sin-0. (19) 

With Eq. @ this yields 

n 2 z => ( 7 ^cff) 4 sin 2 V[l - n 2 z ,(t)}/2, (20) 

and integration of Eq. ( |l"g| ) using n z >(t) = tanh(r) leads 
exactly to Eq. (|l4|). For a ^ 1 this method shows a 
breakdown of Eq. ( |Io| ) because of the impossibility to 
treat the spin-photon interaction as a perturbation. 

In the Landau-Zener experimental setup, the evolution 
of the system described by Eq. (|l(i| ) proceeds for a <C 1 
in two stages, see Fig. |[ The first stage is the LZ pro- 
cess that provides the fraction of molecules P given by 
Eq. (|l3|) in the excited states at t > (the upper energy 
branch in Fig. [|) . In the second stage these excited states 
decay due to the superradiance onto the lower branch in 



0.8- 



0.6 



0.4 



0.2 



0.0 









£ = 0.1, a = 
\e = 0.1, a = 0.01 




a = 0.01 



-40 -20 



20 



40 



60 



80 100 120 
W/A = hvt/A 



FIG. 3: Time dependence of the probability to occupy the 
initial state m during the linear energy sweep with different 
sweep rates and damping coefficients a obtained by numerical 
solution of Eq. (ph|). 



Fig. |l|. According to Eq. (p"3|), at a small sweep rate 
almost all molecules follow the lower energy branch, so 
that the corresponding radiation power is small. To pro- 
duce a powerful coherent electromagnetic radiation the 
resonance should be crossed at a rate v > A 2 /h. To 
obtain the radiation at a fixed frequency, a short longi- 
tudinal field pulse is needed to drive a crystal of molec- 
ular magnets into an excited state with n z i (0) = — cos %j) 
and 71^/(0) = $mip. The time evolution of this state is 
described by Eqs. jT7| ) and (16), whereas the electromag- 
netic radiation is given by Eq. (|l8|). 

Let us now analyze limitations on the superradiance 
relaxation rate of Eq. (|l5|). Firstly, the linear size of 
the sample L should not exceed the wave length of the 
emitted photons A = 27r/fc, otherwise photons emitted 
by different molecules will be out of phase and coherence 
will be destroyed. A conservative estimate is 



kL<l, k = nuo/c = ri A 2 + W 2 /(ch) (21) 
which with N — L 3 /vq yields 



r < Ed A 2 ^ _ {{S z ) o9l i B ) 2 



Ti A 2 + W 2 ' 



(22) 



where Ed is the dipole-dipole energy and vq is the unit- 
cell volume. The maximal rate r max ~ Ed/h and thus 
the maximal intensity of the electromagnetic radiation 
can be achieved at the resonance, W — 0. 

Apart of that, our theory makes use of the weak cou- 
pling a ^ 1 [see comment below Eq. (p0|)]. To work out 
its implications, it is convenient to rewrite Eq. ( |i"4| ) as 



a = (2/3)(fe L) 3 , 



n(A 2 E d ) 1/3 /(ch). (23) 



Thus a ^ 1 converts into k$L < 1. The weak-coupling 
condition is related to the requirement that the dipole- 
dipole interaction is sufficiently weak to prevent the 
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macroscopic system from dipolar ordering at low temper- 
ature (see, e.g., p8| ). In the thermodynamic limit, if the 
condition Ed < A is fulfilled, then in the ground state the 
pseudospins in Eq. (JsJ) point in the x direction instead of 
ordering along the z axis, as favored by the DDI. One can 
see that at the border of the macroscopic region, kL ~ 1 , 
Eq. (||) yields exactly a ~ 1 for E d ~ A and W = 0. 
For Ed < A the weak-coupling condition k$L < 1 is less 
restrictive than the coherence requirement kL % 1 that 
leads to the maximal rate of Eq. (p2]). If A < Ed, then, 
at resonance, the condition k$L < 1 is more restrictive 
than fcL < 1, and the estimation T max ^ Ed/Ti becomes 
invalid. Note, however, that for small tunnel splittings 
A, both k (at resonance) and kg are small and thus both 
conditions are fulfilled for realistic L. 

For Mni2 in the transverse field H x — 6 T, so that 
(S z ) = 0.7775 and A/fi = 2xl0 10 S - 1 (A/k B ~ 
0.15 K > E d /k B ^ 0.04 K), Eq. Q with g = 2, and 
n = I gives a ~ N X 1.9 X 10~ 22 . Taking iV = L 3 /v a 
with the linear size of the crystal L = 0.5 cm and the 
unit-cell volume i>o = 3.7 x 10 -21 cm 3 for Mn.12, one ob- 
tains N = 3.4 x 10 19 and a ~ 0.0063. Then Eq. (§J|) 
gives r ~ 2.5 x 10 s s _1 . For a single crystal of Mni2, the 
maximal possible superradiance rate given by Eq. j2^ ) 
is r max ~7x 10 9 s" 1 . For A/k B ^ 0.15 K this rate is 
achieved at L ~ 1/fc ~ 1.5 cm which is rather large. In 
this case a ~ 0.17, and the weak-coupling approximation 
should still be applicable. For the crystal size L ~ 0.5 cm, 
the maximal rate corresponds to the emission of photons 
of frequency / = cj/(2tt) ~ 10 GHz, i.e., A/k B ~ 0.46 K. 
Similar estimations could be done for Fes which is a bet- 
ter candidate for observing superradiance since in Fes 
the hyperfine interactions that may cause decoherence 
are weak [[15). Taking for Fe 8 v — 2.0 x 10~ 21 cm 3 , 
one obtains r max ~1.3x 10 10 s^ 1 . We have to stress that 
superradiance can be observed for much smaller tunnel 
splittings A than we used in the estimations above, with 
accordingly lower relaxation rates. The only lower lim- 
itations on A are that (i) it has a narrow distribution 
throughout the crystal and (ii) F of Eq. (|l5|) exceeds the 
rates of various incoherent decay processes. 

In conclusion, we have demonstrated that crystals of 
molecular nanomagnets can exhibit superradiance in the 
broad frequency range. Our theory is based upon fun- 
damental electrodynamics of a large spin and is insensi- 
tive to the details of the crystal structure as long as the 
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speqLruin oi molecular spin levels consists oi narrow 
determined by the crystal field and/or external magnetic 
field. Experimentalists should, therefore, focus on large 
single crystals of weakly interacting molecular magnets 
with weak hyperfine interactions. 

Our thinking on this problem has been provoked by 
Javier Tejada. This work has been supported by the 
NSF Grant No. 9978882. 
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